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Force-free magnetosphere of an aligned rotator with 
differential rotation of open magnetic field lines 



O 
O 

(N 

Oh 

<D 

t> 
> 

in 
o 

O 



9 L,: 
6 

CO 



Received: date / Accepted: date 

Abstract Here we briefly report on results of self-consis- 
tent numerical modeling of a differentially rotating force- 
free magnetosphere of an aligned rotator. We show that 
differential rotation of the open field line zone is signifi- 
cant for adjusting of the global structure of the magneto- 
sphere to the current density flowing through the polar 
cap cascades. We argue that for most pulsars station- 
ary cascades in the polar cap can not support stationary 
force-free configurations of the magnetosphere. 
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1 Introduction 

Aligned rotator with the force-free magnetosphere is be- 
ing considered as a first order approximation for the real 
pu lsar magnetosphere since introduction of the model 
bv lGoldreich &: Julian! l)l969t) . Any pulsar model should 
be tested for this simplest case. Recently substantial 



However, in the pulsar magnetosphere current carriers 
are produced mainly by the electron-positron cascades 
and therefore not every curren t density distribution can 
be realized fsee lTimokhinl2006l hereafter Paper I). More- 
over, in Paper I it was shown that the current density 
distribution supporting the magnetospheric configura- 
tion frequently used in theoretical pulsar models, that 
assume the closed field line zone extending up to the 
Light Cylinder (LC) and the open field lines rotating 
with the same angular velocity as the NS, could not be 
realized regardless of a particular model of the polar cap 
cascades. Such models require presence of an electric cur- 
rent flowing against the accelerating electric field in some 
parts of the pulsar polar cap, which cannot be naturally 
explained. 

For the natural stationary configuration of the mag- 
netosphere, when the last closed field line lies in the equa- 
torial plane and magnetic field lines become radial at a 
large distance from the LC (the so-called Y-configuration) , 
the system has two physical "degrees of freedom" , namely 



progress in modeling of pulsar magnetospheres was achieved;) the size of the closed field line zone and ^ the angu _ 
The magnetosphere of an aligned rotator was modeled lar velocity of rotation of the open field lines. In Paper I 
using stationary llContopoulos et. all99alGoodwin et alJ the angu i ar velocity of the magnetic field lines was fixed, 
2004; Gruzinov 20MlTimokhirl2006f) and time-dependent but the size of the closed field line zone was varied. In 
<lKomissarovl200aiBucciantini et al. 2006; McKinney|200a any f the obtained configurations the current density is 
|SpitkovskY||2006D codes. The structure of the magneto- much less than the Goldreich-Julian (GJ) current den- 
sity along field lines passing close to the boundaries of the 
polar cap (see Fig. 5 in Paper I). However, the current 
density, which could be produced by stationary polar cap 
cascades operating in space charge limited flow (SCLF) 
regime, can not be significa ntly less than the correspond - 
ing GJ current density (e.g. lHarding k, Muslimovl El998). 
On the other hand, the current density flowing in the 
polar cap of pulsar having part ially screened p olar gap 
(PSG) cannot be close to zero llGil et al.ll20mh . Hence, 
aligned pulsar with stationary polar cap cascades oper- 
ating in SCLF or PSG regime cannot have stationary 
force-free magnetosphere rotating with a constant veloc- 
ity. As stationary SCLF and PSG are being considered 
as the most probable regimes for the polar cap of pulsar, 



sphere has been obtained even for an inclined rotator 
l)Spitkovskvi2006f) . In these works the angular velocity of 
plasma rotation was assumed to be constant, although 
the case when the open field lines rotate with a con- 
stant angular velocity different from the angular velocity 
of the N eutron Star (NS) was addressed in some work s 
(see e.g. IContopoulosll2005l iBeskin fc Malvshkirjll99^ . 
It was also implicitly assumed that the current density 
in the magnetosphere could adjust to any distribution 
required by the global structure of the magnetosphere. 
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it is worth to find configuration of the force-free magne- 
tosphere compatible with them. 

Here for the first time we consider self-consistently 
the case of a differentially rotating force-free magneto- 
sphere of an aligned rotator in Y-configuration, i.e. we 
are able to explore all possible configurations of such 
magnetosphere. Our aim is to study the impact of dif- 
ferential rotation on the current density distribution in 
the force-free magnetosphere of an aligned rotator. We 
also try to find combinations of the angular velocity dis- 
tribution and the size of the closed field line zone, which 
could be compatible with stationary polar cap cascade 
models. 



2 Equations for differentially rotating force-free 
magnetosphere 

Magnetic field lines in the closed zone are equipotential 
and plasma there corotates with the NS. In the open 
field line domain the angular velocity of plasma Sip is 
different from the angular velocity of the NS rotation fl. 
The angular velocity of the open magnetic field lines Qp 
is determined by the potential difference in the polar cap 
of pulsar. In stationary cascade models flp is close to Q 
for young pulsars (see Sec. 2.2 in Paper I.) 

The force-free magnetosphere of an aligned rotator 
with differential rotation can be described by a solu- 
tion o f the so-called pulsar equation, derived bv lokamotol 
l)l974t) . In notations of Paper I (see eq. 20 there) it has 
the form 



(3 2 x 2 



I 



_ s d_s 



d x ip- 



x ^^± (vvo =o. 



dip ' dip 



(1) 



Here /3 = fip /£2, S and ip are normalized poloidal current 
and magnetic flux functions correspondingly, z axis is co- 
aligned with fi. All coordinates are normalized to the LC 
radius of a corotating magnetosphere i?LC = c/Q. The 
last term in the equation takes into account contribu- 
tion of differential rotation into the force balance across 
magnetic field lines. The current density in the open field 
line zone of the magnetosphere depends on differential 
rotation through condition at the true LC, which is at 
x(ip) = c/i?p(V0 (see Sec. 2.3 in Paper I) 



(2) 



Changes of (3 results in changes of the poloidal current 
density j being proportional to dS/dip. 

So, the differential rotation contributes to the force 
balance perpendicular to the magnetic field lines, changes 
the position of the LC and affects the current density dis- 
tribution, required for smooth transition of the solution 
trough the LC. Although, the differential rotation can 
modify the current density distribution, it is a quantita- 
tive question whether or not a given differential rotation 




Fig. 1 Structure of the magnetosphere for xo = 0.7 and the 
angular velocity of the open magnetic field lines Qp{ip) = 
J?/3(i/>) (see text). Magnetic field lines are shown by the thin 
solid lines, the last closed field line - by the thick solid line, 
the position of the Light Cylinder is shown by the dot-dashed 
line. 



f2p{ip) and the corresponding current density distribu- 
tion j(tp) could be supported by stationary cascades in 
the polar cap of pulsar. We note that f3(ip) is a free pa- 
rameter in this problem, while j(tp) is obtained from the 
solution. 

For solution of equation we developed an ad- 
vanced version of the multigrid code described in Pa- 
per I. Now the position of the Light Cylinder is found 
self-consistently on each iteration step and poloidal cur- 
rent is calculated from the equation at the LC. The 
return current flowing in the current sheet between the 
closed and open field line zones is smeared over the inter- 
val ["0iast, V'last +dip] 1 , and the angular velocity is contin- 
uously changing in the same interval from some value at 
the last open filed line ^(V'last) to fl at ip = "0iast + dip. 
Equation Q is solved in the whole domain including the 
current sheet. Smearing of the return current and con- 
tinuous changing of fip in the current sheet allows us to 
take into account the contribution of the current sheet 
into the force balance across magnetic field lines at the 
boundary between closed and open field line zones. 



3 Main results 

Stationary polar cap cascades without particle inflow 
from the magnetosphere would provide more or less con- 
st antcuT£ejrt_dCTsitv_djste cap (see 
e.g. lHarding fc Muslimovl Il998t ICil et al] 12003^ . How- 
ever, in configurations of the force-free magnetosphere 
with constant fip the current density goes to zero near 

1 '(/'last is the value of the normalized magnetic flux function 
corresponding to the field line separating the closed and open 
field line domains 
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Fig. 2 Current density distributions in the polar cap of pul- 
sar jp C (d) corresponding to ~ j (see text) are plotted 
for different xq as functions of the colatitude 8. j pc is nor- 
malized to the Goldreich- Julian current density \jaj\. 8 is 
normalized to the colatitude of the polar cap boundary 8 pc . 
The current d ensity corresponding to the Michel's solution 
(Michel 1973) is shown by the dashed line. The distribution 
for xo = 0.1 practically coincides with the distribution for 
the split-monopole case. 



the polar cap boundaries. Here we try to find configura- 
tions of a differentially rotating force- free magnetosphere 
with current density, which is nearly constant over the 
polar cap at the NS surface. Current density in the mag- 
netosphere rotating with a constant angular velocity de- 
creases toward the polar cap boundary because of the 
boundary condition ip(x > xq) — ipust {xq is the po- 
sition of the point, where the last closed field line in- 
tersects the equatorial plane). In the equation for the 
poloidal current (eq. J2J) the first term on the left hand 
side goes to zero for ip approaching V'last- Hence, in or- 
der to increase the current density near the polar cap 
boundaries d(3/dip must be positive, i.e. the maximum 
value of p should be achieved at the polar cap boundary, 
for ip — ^iast' Obviously, the maximum possible value of 
j3 is 1. Such behavior of (3 is expected in models of the 
polar cap cascades too. 

We have performed computations for different sizes 
of the closed field line zone x$. In each of these simula- 
tions the distribution of (3 over ip was adjusted ad hoc 
in order to obtain nearly constant current density over 
the polar cap, which approaches zero only at field lines 
very close to the polar cap boundary 2 . For any size of 
the closed field line zone Xq it is possible to construct a 
set of force-free magnetospheric configurations with the 
current density distribution being almost constant over 
the polar cap, j(ip) — j = const by choosing different 

One of the important properties of the obtained set of 
solutions is that for each xq the constant current density 
j cannot be made greater than some maximum value j. 

2 we assume that the current sheet carries only the return 
current, so j changes the sign at V'last- 




Fig. 3 Angular velocities of rotation of the open field lines 
O f (8) — Q /3(8) (see text) for different xq as functions of the 
colatitude 8 in the polar cap of pulsar. 8 is normalized to the 
colatitude of the polar cap boundary 8 pc . Qf{8) is normalized 
to the angular velocity of the NS rotation Q. The distribution 
for xo = 0.1 practically coincides with the distribution for the 
split-monopole case. 



Distribution j(ip) ~ J corresponds to a distribution /3(tp), 
which achieves 1 at the polar cap boundary, /3(V>iast) = 1. 
As it was stressed above, the maximum value of (3 is 
achieved at -0 = ^last- So, (3 differs from other distri- 
butions /3(V0 which provide constant current density, in 
that it achieves the maximum possible value. Both j and 
/3(V0 depend on xq. 

The differential rotation in the magnetosphere of a 
young pulsar with stationary polar cap cascades is rather 
small (see Section. 2.2 in Paper I). So, solutions with 
the minimal deviation of (3(ip) from 1 would be of most 
interest for us. These turned out to be solutions with 
(3 = /3(V0i i- e - with j ~ j. One of such magnetospheric 
configurations with xq — 0.7 is shown in Fig. ^ Impor- 
tant properties of solutions with j — j are: 

— The value J increases with decreasing of xq , however 
it does not exceeds Jgj/2, the half of the GJ current 
density near the NS surface (see Fig.EJ). 

— The angular velocity of rotation of the open mag- 
netic field lines has smaller deviation from fl with 
decreasing of xo, see Fig.|3] As in the case of the cur- 
rent density distribution there exist an asymptotic 
distribution of fi-p(ip). This asymptotic distribution 
deviates strongly from ft(ip) = 1. 

— The asymptotic form of /3(tp) and j for xq — > matches 
the corresponding distributions in the split monopolc 
configuration. 

— Energy losses of the aligned rotator increases with 
decreasing of xq. This increase goes a bit faster than 
in the case of constant fl-p , cf . Fig. 0] here and Fig. 8 
in Paper I. 

— The total energy of the magnetosphere increases with 
decreasing of xq. 

— In none of the considered configurations the force- 
free condition is violated in the calculation domain, 
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cap cascades. In this case the characteristic current den- 
sity flowing through the gap must be also much less than 
the GJ current density. 
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Fig. 4 Energy losses of aligned rotator normalized to the 
magnetodipolar energy losses as a function of xo for differ- 
entially rotating magnetosphere with J2f(V0 = ^/^WO ( see 
text). 



i.e. the electric field is everywhere smaller than the 
magnetic field. 

Let us now discuss implications of these solutions for 
the physics of pulsars. 



4 Discussion 



If the polar cap cascades operate in stationary SCLF 
regime the current density in the polar cap of pulsar 
near the NS surface is very close to the local GJ current 
density. However, due to inertial frame dragging (i.e. due 
to changing of the effective J?) the local GJ charge den- 
sity decreases slower with the distance than the charge 
density carrie d by the charge-separated flow from the 
NS surface fsee lMuslimov fc Tsvganll99l . This discrep- 
ancy achieves ~ 15% at a distance of several radii of the 
NS. So, the current density at that distance will be by 
~ 15% less then the corresponding local GJ current den- 
sity. Here we solve the problem in flat space-time. The 
boundary conditions on the NS surface and, hence, the 
physical quantities we used for normalization should be 
taken at some distance from the NS, where GR effects 
are negligible. The GR corrections will result in ~ 15% 
higher ratio of j to the local value of j G3 at the NS sur- 
face than the one obtained in our solution. However as 
it follows from our results, there is no configurations of 
the magnetosphere with almost constant current density 
larger than j G3 /2, Hence, stationary force-free magne- 
tosphere of an aligned rotator can not be supported by 

stationary polar cap cascades oper ating in SCLF regime. , , 

Po lar cap cascades operating in lRuderman &, SutherlanflhiTnokhii] l)200 
l|l975l) regime with almost vacuum electric field in the 
gap allow larger deviations of the current density from 
j GJ as well as larger deviation of (3 from 1. For older 
pulsar, when (3 could be rather small, even configuration 
with xq approaching 1 could be supported by the polar 



F or the partially screened polar gap (PSG) l(Gil et alJ 
2003) the presence of a strong multipolar component of 
the magnetic field in the polar cap of pulsar is essen- 
tial. In this model the current density at the NS surface 
for young pulsars should be close to the local value of 
j GJ . However, the local GJ charge density is determined 
by the non-dipolar magnetic field. For this field the an- 
gle between B and fl can be rather large and the GJ 
charge density will be accordingly less. At some distance, 
where magnetic field becomes dipolar and the angle be- 
tween B and fl is small, the current density of the flow 
would become smaller than the local GJ current density. 
This effect can be much stronger than the difference be- 
tween j and j GJ due to inertial frame dragging. Hence, for 
some specific configurations of the local magnetic field 
at the NS surface, stationary polar cap cascades could 
support a stationary force-free configuration of the mag- 
netosphere. 

So, aligned pulsars with polar cap cascades operat- 
ing in SCLF regime, young pulsars with "Ruderman- 
Sutherland" cascades as well as pulsars with PSG and ar- 
bitrary surface magnetic field configuration can not have 
both polar cap cascades operating in stationary regime 
and a stationary force-free magnetosphere. From ener- 
getic point of view, configurations with non-stationary 
polar cap cascades would seem to be preferable. The en- 
ergy of the magnetosphere increases with decreasing of 
Xo- The magnetosphere would try to achieve the config- 
uration with the smallest possible energy, i.e. with the 
largest possible xo, compatible with the physical con- 
ditions set by the polar cap cascades. If there are sev- 
eral possible self-consistent configuration of the magne- 
tosphere and polar cap acceleration zone, the configura- 
tion with the largest xowill be realized. With aging of 
the pulsar the conditions in the cascades change. This 
would result in changing of the magnetospheric configu- 
ration and, hence, in changes of energy losses relative to 
the corresponding magnetodipolar losses. This yields the 
breaking index different from 3. However, more detailed 
study of polar cap cascade properties is necessary in or- 
der to construct self-consistent magnetospheric configu- 
rations and to obtain the value of the breaking indexes. 
The current adjustment is necessary also for inclined ro- 
tator and in this case similar effects would be present 
too. 



More detaile d description of the results is given in 
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